EUROPHYSICS LETTERS 
Europhys. Lett., (), pp. () 



Higher moments of spin-spin correlation functions for the 
ferromagnetic random bond Potts model 

Marc- Andre Lewis (*) 

Laboratoire de Physique Theorique et des Hautes Energies, (**) 

Universites Pierre et Marie Curie (Paris VI) et Denis Diderot (Paris VII). 

Boite 126, Tour 16, ler etage, 

4 pi. Jussieu, 75251 Pans CEDEX 05, FRANCE 

(received ; accepted ) 

PACS. 64.60Ak - Renormalization group studies of phase transitions. 
PACS. 11.25Mj - Conformal field theory; algebraic structures. 



Abstract. - Using conformal field theory techniques, we compute the disorder-averaged pth 
power of the spin-spin correlation function ((a(0)a(R)) p , p 6 Z) for the ferromagnetic random 
bond Potts model. We thus generalize the calculations of Dotsenko, Dotsenko and Picco, 
where the case p = 2 was considered, and of Ludwig, where first-order computations where 
made for general p. Perturbative calculations are made up to the second order in e (e being 
proportional to the central charge deviation of the pure model from the Ising model value). The 
explicit dependence of the correlation function on p gives an upper bound for the validity of the 
e-expansion, which seems to be valid, in the three-states case, only if p < 4. 



Since the first calculations made by Ludwig || , a lot of attention was given to the study of the 
random bond Potts model. It was established that the introduction of randomness changes 
the critical behaviour of the system, as predicted by the Harris criterion. Using perturbative 
conformal field theory techniques Q and e-regularization which consists here in a shift of 
central charge from the Ising model value, first order || and second order calculations [|| 
clearly established the existence of fixed points in the renormalization group flow. In fact, 
there exists two different fixed point solutions: one with replica symetry (RS) and another 
where this symetry is broken (RSB). Recent results by Dotsenko, Dotsenko and Picco [[[] 
support the RS fixed point critical behaviour of the random bonds Potts model. To compare 
both schemes, they compute the disorder-averaged second moment of the spin-spin correlation 
function (a(0)a(R)) 2 with broken and unbroken replica symetry. Numerical simulations for 
the 3-states and 4-states models don't shown significant deviation from the replica symetric 
solution. 
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However, the observed higher moments of correlation functions seem to be in contradiction 
with the values predicted by the RS e-expansion calculations [Q. In this letter, we will compute 
the disorder-averaged p-th power of the spin-spin correlation function (((t(0)<t(R))p) in the 
replica symetric case. The explicit dependence of this quantity on p shows how the expansion 
validity breaks down for sufficiently large p. We find, for the 3-state Potts model, that the 
expansion is valid only if p < 4, thus confirming the difference between observed and predicted 
values for high moments. 

The partition function of the nearly-critical g-states random bond Potts model, is well 
known to be of the form 

Z{0) = Tr exp{-/f - (1) 

where Hq is the Hamiltonian of the conformal field theory corresponding to the g-states Potts 
model with coupling constant Jo the same for each bond. The Hamiltonian H%, being the 
deviation from the critical point induced by disorder is of the form 

Hi = J d 2 xr{x)e(x), (2) 

where t(x) ~ f3J{x) — (3 c Jq is the random temperature parameter. The theory is defined on 
the whole plane. We shall assume, for simplicity, that t(x) has a gaussian distribution for 
each x, with 

P~Pc 



t(x) = t = — — (3) 

(t{x) - t q ){t(x>) - t q ) = g 6W(x-x') (4) 

The usual way of averaging over disorder is to introduce replicas, that is, n identical copies 
of the same model for which: 

n r, n 

{Z{P)) n = Tr exp{-]T/4 a) - / Ar(^ £a (x)}. (5) 

a=l J a=l 

Taking the average over disorder, by performing gaussian integration, one gets 

n r, n r, n 

(Z(P)r = Tr exp{-^/4 a) -T / d 2 x £ e a (x) + g / d 2 x £ e a {x)e b (x)}. (6) 

a=l J a=l a^b 

This is a field theory of n coupled models with coupling action given by 

„ n 

H \xit = / d 2 x^2e a (x)e b (x). (7) 

Note that only non-diagonal terms are kept since diagonal ones can be included in the Hamil- 
tonian Hq. Moreover, they can be shown to have irrelevant contributions, since their OPE 
consist of the identity plus terms that are irrelevant at the pure fixed point. We now turn our 
attention to the p-th moment of the spin-spin correlation function (a(0)a(R)) p . In terms of 
replicas, it can be written as 

= Um^-f^ Y, M°)"--M0) E v bl (R)---v bp (R))- (8) 
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The operator to be renormalized is then 

O p (x) = a ai (x)a a2 (x) ■ ■ ■ a ap (x), ai ^ a 2 ■ ■ ■ ^ a p , 1 < a t < n, (9) 
perturbed by the interaction term; 

O p (x) = O p exp{-H int } = O p (l- H- mt + \{H int ) 2 (10) 

We will define the amplitude Z, for which we will derive RG equations, as 

6 p (x) = ZO p (x). (11) 

The task at hand is thus to rewrite O p (x) in the form ([ll]) by doing all possible contractions 
and operator algebra. We will compute Z up to the second order in g . To do so, we use 
the Coulomb gas formulation of minimal conformal held theories Q. In this formalism, the 
central charge of the theory is written as 

c = 1 - 24al (12) 



a(x)e(y) = — r^-a(x) + finite contributions, (15) 



a± = a ±ya 2 + l aq_a_ = — 1. (13) 

For the Ising model, a 2 + — 4/3 and c = 1/2, while for the 3-states Potts model, a 2 + — | and 
c = |. For a generic model, we will write a\ = | — e, so that c = i + ^-e + C(e 2 ). In particular, 
e = corresponds to the 3-states Potts model. For Potts models, the energy operator e{x) is 
the primary field <&x,2 so that its conformal dimension is 

A, = A ll2 + A li2 = (a - +2Q + )2 - (Q - +a+)2 = 1 - \e. (14) 

We shall often use the spin and energy operators product expansion 

D 

\x - y\ l 

where D, the operator algebra coefficient, is known to be \ + 0(e 2 ) @]. One can get rid of the 
finite terms by projecting correlations functions on <r(oo). 

Renormalization group equations will be derived by integrating from a cut-off of 1 (in lattice 
spacing units) to a new one a (a 3> 1). First order calculations are straightforward. Since 
operators with different replica indexes have zero product expansion, there are p(p— 1) possible 
contractions, that is 

/n 

-> a ai (x)---a ap (x)gp(p-l) d 2 y{a(x)e(y)a(oo)) 2 

Jl<\y~x\<a 

f D 2 

- a ai (x)---a ap (x)gp(p-l) dy 

Jl<\y-x\<a \ x ~ V\ 

= O p (x)p(p-l)g^a 3 \ (16) 



3e 

So, the first order corrections to Z are 



5ZU = Zp(p-l)g^a 3 *. (17) 
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Second order calculations require more work. There are five different types of contractions 
possible; three of them occur if p > 2, the fourth if p > 3 and finally the fifth if p > 4. The 
first three diagrams were computed in Q, and, for generic p, only combinatorial factors are 
modified. We will only give their expression and concentrate on the computation of the two 
last diagrams. The first diagrams give the contributions 

D {2) = O p (x)p(p-l)(n-2)g 2 ^^(l + eK)a^ (18) 

(1) 

D% } = finite contributions (19) 
Df = 0p (x)p(p-ltf(^-£)a* (20) 

where K — 6 log 2. We only consider the divergent part of the diagrams since these are the 
only ones appearing in the RG equations. The fourth diagram expression is given by 

= O p (x)j-^g 2 J |d 2 yd 2 y'(a(0) £ (y) e ( 2/ ')^(oo)}(a(0) e ( 2/ )a(oo))(a(0)e( 2/ ')^(^)) 
= O p (x)j-^-g 2 D 2 J J d 2 yd 2 y'\y\^\y'\-^{a{Q)e(y)e{y') < j(^)). (21) 

A trivial change of variable and the use of the fact that 

(a(0)e(y)s(y')a(oo)) = A 2Ae (a(0)e(Xy)e(Xy')a(oo)) 

leads to 

D {2) = 20 p (x)p(p~l)(p-2)g 2 f d 2 y'\y'\ 2 -^ [ d 2 y\y\- A ' (*(0)e(l)s(y)*(oc)) 



= 2nD 2 — J d 2 y\y\-^(a(0)e(l)s(y)a(^)). (22) 

The calculation of this integral is done with the use of the techniques described in |jj . One 
gets 

= O p (x)p(p-l)(p~2)g 2 ^(8 + ea)a 6 \ (23) 

with a = 33-22^. 

(2) 

The calculation of is simpler. The diagram consists of four ae contractions: 
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4 2) = \o p {x)p(j>-l){p-2)(j>-Z)g 2 ^J d 2 y^ 

= O p (x)p(p-l)(p-2)(p-3)g 2 ^^a^ (24) 
Collecting all results, we get the second order correction to Z: 



/ a 2 r>2 

6Z& = Zg 2 p(p-l)a 6e Un- 2)^-^(1 + eK) + 



47r^ 4 

Ik 2 36e 



'it 2 D 2 2tt 2 D 4 
^(P-2)(^-(8 + ea) + (p-3)-^-) ) (25) 
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We can now write the RG equation for Z (£ = logo): 



where 



J = a^- = Z (A(p, e)g(a)a 3 * + B(p, e)g 2 (a)a 6e ) , (26) 



A{p) = 2irD 2 p(p-l) (27) 

bm - P& >- d((»- 2)?^(i + do + (5^-: 

+p(p - l)(p - 2) f —(8 + 6a) + (p - 3)^^j ■ (28) 

There is also a renormalization of the coupling constant g. Calculations were originally 
presented in j^]; we shall not review them here. For a given cutoff a, g renormalizes as (tilded 
operators will represent renormalizcd quantities) 

.g = a 3e (.9 + W^), (29) 

with the cut-off dependent factor introduced to return to the cut-off scale a = 1 . We now 
invert the renormalization equation up to the second order in g: 

9 = a- 3 T9-^f) (30) 

~ f AD 2 \ 
Z=Z\l-—p(p-l)g\. (31) 

Replacing bare quantities by renormalized ones in (p6~|), and using the fact that D = ~ + 0(e 2 ), 
one gets (we let g 

d ~W = ^ )pip - 1} (l m + ( (n - 2) ^ K - k + ^' 2) W2 a ) ^) < (32) 

where, we recall, K = 6 log 2 and a = 33 - 

We can now easily solve the RG equation (p2[). ft can be rewritten in the form (dropping 
the tildes), 

^p- = liOHO (33) 

7(0 = p(p-l)(~g(0+((n-2)±K-± + (p-2)±a^g 2 (0Y (34) 

To compute the correlation functions, it will be useful to assume the RG evolution to go 
from the lattice cut-off (~ 1) to the scale R (we write £r = log R). To do so, we need the fixed 
point value of g, which we will note g*. It is know to be of the form J3|, |^]: 

^ = ^+ 9 r 2 + 0(e 3 ). (35) 

Taking the limit on the number of replicas (n — 0) and using the explicit form of one 
obtains the fixed point value of 7, noted 7* 

* = >-'>(Hs-x + sM'°) +0 <«'>- < 36 > 
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We are now able to compute the correlation functions. Using scaling laws, we get 
(a(0)a(R))P = lim^fq £ a ai (0) • • • a ap (0) £ a bl (R) ■ ■ ■ <r bp (R)) 



lim 



(^(^)) 2 



(37) 



The final result is obtained by using the fixed point value Z(£r) ~ e 7 *^ H = R 1 * . One thus 
gets 



1 

with 



a(0)a(R))P ~ T ^ 7 -. (38) 



A' ffP =pA ff - 7 *. (39) 

The deviation from the pure model is thus given by 7*. Having show this quantity to be 
of the form Ae + Be 2 + 0(e 3 ), we can now look at the domain of validity of the e-expansion. 
Evidently, it becomes absurd if |Ae+i?e 2 | ~ |pA CT p |. For the 3-states Potts model, this happens 
for p > 5. This explains why this method cannot predict disorder- averaged moment for such 
p's. In contrast the expansion makes a good approximation for p < 4. 

To conclude, let us derive another physically interesting quantity, which is the derivative 
of A^ P with respect to p, evaluated at p = (it is ao (not to be confused with the Coulomb 
gas parameter) in Ludwig's notation). It describes the asymptotic decay of the spin-spin 
correlation function ((a(0)a(R)) ~ R >} c , ) . It is straightforwardly shown to be 



dA> 




This quantity is probably the easiest to measure in numerical simulations. 
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